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Given an undirected graph with edge costs and a subset of k > 3 nodes called terminals, a multiway, 

(N 

J> ' or fc-way, cut is a subset of the edges whose removal disconnects each terminal from the others. The 

multiway cut problem is to find a minimum-cost multiway cut. This problem is Max-SNP hard. Recently 
Calinescu, Karloff, and Rabani (STOC'98) gave a novel geometric relaxation of the problem and a 

(N ' 

rounding scheme that produced a (3/2 — 1/ ^-approximation algorithm. 

In this paper, we study their geometric relaxation. In particular, we study the worst-case ratio between 
the value of the relaxation and the value of the minimum multicut (the so-called integrality gap of the 
relaxation). For k = 3, we show the integrality gap is 12/11, giving tight upper and lower bounds. That 
is, we exhibit a family of graphs with integrality gaps arbitrarily close to 12/11 and give an algorithm 
that finds a cut of value 12/11 times the relaxation value. Our lower bound shows that this is the best 
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possible performance guarantee for any algorithm based purely on the value of the relaxation. Our upper 
bound meets the lower bound and improves the factor of 7/6 shown by Calinescu et al. 

For all k, we show that there exists a rounding scheme with performance ratio equal to the integrality 
gap, and we give explicit constructions of polynomial-time rounding schemes that lead to improved 
upper bounds. For k = 4 and 5, our best upper bounds are based on computer-constructed rounding 
schemes (with computer proofs of correctness). For general k we give an algorithm with performance 
ratio 1.3438 - e k . 

Our results were discovered with the help of computational experiments that we also describe here. 

1 Introduction 

As the field of approximation algorithms matures, methodologies are emerging that apply broadly to many 
NP-hard optimization problems. One such approach (c.f. HJ |9] E El ED) has been the use of metric and 
geometric embeddings in addressing graph optimization problems. Faced with a discrete graph optimization 
problem, one formulates a relaxation that maps each graph node into a metric or geometric space, which in 
turn induces lengths on the graph's edges. One solves this relaxation optimally and then derives from the 
relaxed solution a near-optimal solution to the original problem. 

This approach has been applied successfully [ 3 1 to the min-cost multiway cut problem, a natural gener- 
alization of the minimum (s,i) -cut problem to more than two terminals. An instance consists of a graph 
with edge-costs and a set of distinguished nodes (the terminals). The goal is to find a minimum-cost set of 
edges whose removal separates the terminals. If the number of terminals is k, we call such a set of edges a 
k-way cut. 

The first approximation algorithm for the multiway cut problem in general graphs was given by Dahlhous, 
Johnson, Papadimitriou, Seymour, and Yannakakis |5). It used a traditional minimum (s,t)-cut algorithm 
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as a subroutine and had a performance guarantee of 2 — 2/ k. 

In the work that prompted ours, Calinescu, Karloff, and Rabani used a novel geometric relaxation 
of fe-way cut in a (3/2 — l/£:)-approximation algorithm. Their relaxation uses the k-simplex A = {x € 
R fc : x > 0, J2i x i = 1}> which has k vertices; the i th vertex is the point a; in A with X{ = 1 and all other 
coordinates 0. The relaxation is as follows: map the nodes of the graph to points in A such that terminal i is 
mapped to the i th vertex of A. Each edge is mapped to the straight line between its endpoints. The goal is 
to minimize the volume of G, 

vol(G) = ^2 cost ( e ) • M 

edges e 

where cost(e) is taken to be the cross-sectional area of edge e and \e\ denotes the length of the embedded 
edge e, defined as half the L\ distance between its endpoints. The problem of finding an embedding that 
minimizes the volume can be formulated as a linear program (LP). The factor half in the length function is 
present to scale the distance between terminals to 1, so the LP is a relaxation of the minimum /e-way cut 
problem. 

To see that the given LP is a relaxation of fc-way cut, consider any /c-way cut and let Si be the set of 
nodes reachable from terminal i in the graph with the cut-edges removed. Consider a geometric embedding 
in which all nodes in Si are mapped to vertex i of A. For any edge, its embedded length is either 0, if the 
endpoints lie in the same Si, or 1, if the endpoints lie in distinct Si. Hence the volume of this embedding is 
equal to the cost of the fe-way cut. 

The algorithm of Calinescu et al. finds a minimum volume embedding by linear programming. It then 
uses a randomized rounding scheme to extract a cut from this embedding. Ignoring the graph, the scheme 
chooses (from a carefully selected distribution) a k-way cut of the simplex — a partition of the simplex into 
k polytopes, each containing exactly one vertex of the simplex. The fe-way cut of the simplex naturally 
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induces a fc-way cut in the embedded graph — namely, the set of edges with endpoints in different blocks of 
the partition. This cut has expected cost at most 3/2 — 1/ k times the volume of the embedding. 

1.1 Our results 

Our goal is to further understand the geometric relaxation, with the hope of developing better approximation 
algorithms. We aim to determine the integrality gap of the relaxation and to find an algorithm whose 
approximation ratio matches the integrality gap. (Formally, the integrality gap is the supremum, over all 
weighted graphs G, of the minimum cost of any /c-way cut of G divided by the minimum volume of any 
embedding of G. Note that the integrality gap is the best approximation ratio we can prove using an analysis 
that bounds the optimum cut only by the value of the relaxation.) 

In this paper, we resolve this question for 3-cut and provide improved results for the general k-cut 
problem. For k = 3 we give a rounding algorithm with performance ratio 12/11, improving the Calinescu 
et al. bound of 3/2 — 1/3 = 7/6. We also show that 12/11 is the best possible bound, exhibiting a graph 
family with a gap of 12/11 — o(l) between its embedded volume and minimum 3-way cut. Thus, for k = 3, 
we determine the exact integrality gap and give an optimal rounding algorithm. 

For larger k, we obtain results based on both computation and analysis. We give a non-constructive proof 
that, for every k, there exists a (not necessarily polynomial time) rounding scheme whose performance 
guarantee equals the integrality gap. For k = 4,5, we use LP-derived and -analyzed rounding schemes 
to give explicit approximation bounds of 1.1539 and 1.2161 respectively, improving the corresponding 
Calinescu et al. bounds of 1.25 and 1.3. For larger k we give an algorithm obtaining a (analytic) bound of 
1.3438 — €k where e k > 0. The quantity e k can be evaluated computationally for any fixed k; we use this to 
prove that 1.3438 - e k < 3/2 - 1/k for all k. 

Our efforts to find geometric cutting schemes that achieve good guarantees were guided by experiments: 



4 



we formulated the problem of determining an optimal probability distribution on /c-way cuts of the simplex 
as an infinite-dimensional linear program and solved discrete approximations of this linear program and its 
dual. From these solutions we were able to deduce the lower bound and, using that, the upper bound for 
k = 3. These experiments also guided our search for cutting schemes that work for larger values of k. 
The upper and lower bounds for k = 3 were discovered independently by Cunningham and Tang [4 1. 

Presentation overview. In Section|2]we discuss the geometric ideas underlying the problem. In Section|3] 

we describe the computational experiments we undertook and the results it gave for small k. In Sections |4] 
and |5] we resolve the 3-terminal case, giving matching upper and lower bounds. Finally, in Section |6] we 
present our improved algorithm for general k. In the appendix we prove that, for all k, there exists a rounding 
scheme matching the integrality gap. 

2 The geometric problem 

Finding the integrality gap of and a rounding scheme for the relaxation turns out to be expressible as a 
geometric question. That is, we can express integrality gaps and algorithmic performance purely in terms of 
the simplex, without considering particular graphs or embeddings. 

Consider an edge e, which under the relaxation is embedded as a line segment in the simplex. We 
overload e to denote this embedded segment as well. For any segment (or edge) e, we let e£ denote the 
projection of e onto the £ th coordinate axis, namely the one dimensional interval {xi \ x € e}. We write 
min e£ = min xge xg for the minimum value in the projected interval, max = max xge X£ for the maximum 
value in the projected interval, and \eg\ = maxeg — mine^. Finally, as mentioned in Section^ the length 
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|e| of an edge e is defined to be half its L\ norm, that is, 

k 

N=EN/ 2 - 

2.1 Density 

Recall that a fc-way cut of the simplex is a partition of the simplex into k polytopes, each containing a unique 
vertex of the simplex, and that such a cut induces a fc-way cut of any embedded graph. By a cutting scheme, 
we mean a probability distribution P on /c-way cuts of the simplex. For any line segment e, the density of P 
on segment e, denoted Tfc(P, e), is the expected number of times a random cut from P cuts e, divided by the 
length |e| of e. 1 

Define the maximum density of P, Tk{P) and the minimal maximum density rjt as follows: 

Tk (P) = sup Tfc (P, e) and t£ = inf T k (P) , 

e r 

There is always a line segment of infinitesimal length that achieves the maximum density, since any segment 
can be divided into two edges, one of which has density no less than the original. Thus, in the remainder of 
this paper, we will focus discussion on such infinitesimal segments. 

The relevance of t£ is the following (this is implicit in the work of Calinescu et al.): 

Lemma 2.1 For any cutting scheme P and embedded graph G, the expected cost of the k-way cut of G 
induced by a random k-way cut from P is at most Tk{P) times the cost of the embedding of G. 

Corollary 2.2 Any cutting scheme P yields a randomized approximation algorithm with approximation 
ratio at most Tfe(P). 

'Note that in principle a line segment could be cut more than once by the fc-way cut of the simplex. We therefore speak of the 
expected number of times that e is cut, rather than the probability that e is cut. 
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Proof: The endpoints of any edge e are embedded at two points in the simplex, so the edge corresponds to a 
segment connecting those two points. The expected number of times the edge is cut is at most Tfc(P, e) • |ej. 
By the Markov inequality this upper bounds the probability that the edge is cut. Thus, the expected cost of 
the /c-way cut is at most 

2(7*(P,e)-|e|)cost(e) < r k {P) ^ |e| ■ cost(e) 

e 

= t*(P)vo1(G). 

We have already argued that vol(G) lower bounds the optimum fc-way cut, so the result follows. rj 

The above argument implies that no cutting scheme P can have a maximum density Tfc(P) below the 
integrality gap. In fact, we show that there always exists a cutting scheme whose maximum density equals 
the integrality gap. 

Theorem 2.3 There exists a cutting scheme whose maximum density equals the integrality gap, thus, t% 
equals the integrality gap. 

We give the proof in the appendix. The proof is based on the observation that the problem of choosing a 
rounding scheme to minimize the performance ratio is itself a (infinite dimensional) linear programming 
problem; furthermore its dual is the problem of choosing a weighted graph to maximize the integrality 
gap. This observation seems to hold in a fairly general setting beyond the /c-cut problem (details are in the 
appendix). 

Calinescu et al.'s algorithm gives a cutting scheme showing that r£ < 3/2 — 1/k. In this paper we show 
that r 3 * = 12/11, and that, for all k, r* k < 1.3438. 
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2.2 Alignment 

We have just argued that the key question to study is the maximum density of (infinitesimal) line segments 
relative to a cutting scheme. Calinescu et al. showed that one can restrict attention to segments in certain 
orientations. We say a segment e in A is i,j -aligned if e is parallel to the edge connecting vertices i and j 
of A. We say it is aligned if it is i, j-aligned for some pair of vertices. Calinescu et al. observed that the 
endpoints of any segment e can be connected by a piecewise linear path of total length |e| whose segments 
are aligned. The segment e is cut if and only if some edge on this path is cut. Given any embedding of a 
graph, Calinescu et al. apply this transformation to each segment connecting two embedded vertices, without 
changing the volume of the embedding. Thus, they show that without loss of generality one may restrict 
attention to embeddings in which all edges are aligned. 

Fact 2.4 Segment e is i,j-aligned if and only if\e\ = |ej| = \ej\ and \ee\ = Ofor I ^ 
(Note that |e| denotes half the L\ norm, while \e\ \ and |e2 1 are standard absolute values.) 

2.3 Side parallel cuts (SPARCS) 

In this paper, we mainly restrict attention to a particular class of cutting schemes. Define A Xi=p = {x € 
A : Xi = p) and A Xj > p = {x G A : Xi > p). Note that A Xi=p is a hyperplane that runs parallel to the face 
of the simplex opposite terminal i and is at distance p from that face; it divides the simplex into two parts, of 
which A Xj > p is the "corner" containing terminal i. An i, j-aligned segment (x, y) is cut by the hyperplane 
A Xe=p if and only if t G {i, j} and p is between X£ and y^. 
We define a side-parallel cut (spare) of the simplex: 

1. Choose a permutation a of the vertices; 

2. For each vertex i in order by a (except possibly the last), choose some pi G [0, 1]; 
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3. Assign to vertex i all points of A Xi > Pt not already assigned to a previous terminal. We say that 
terminal i captures all these points, and that terminal i cuts an edge e if it captures some but not all of 
the yet-uncaptured part of e. 

This scheme cuts up the simplex using hyperplanes A Xi=p . In this context, we call each A x . =p a slice. 

We consider algorithms that sample randomly from some probability distribution over spares. Our 
restriction to spares was motivated by several factors. The rounding algorithm of Calinescu et al. uses 
only spares. Furthermore, our computational study of the 3-terminal problem (discussed below) and some 
related analytic work gave some evidence that the optimal algorithm was a distribution over spares (this 
conjecture was confirmed analytically for the 3-terminal case). Lastly, spares have concise descriptions (as 
sequences of k — 1 slicing distances) that made them easy to work with computationally and analytically. It 
is conceivable, though, that one might do better with cuts that are not spares. For example, one might wish 
to slice off two terminals simultaneously, and then separate the two from each other. Indeed, we know of no 
proof that for k > 3 the optimal cut must be made up of hyperplanes; curved surfaces might do better. 

For segment e, recall that et is the interval {xi \ x G e}. The key properties of spares are expressed in 
the following fact. 

Fact 2.5 An i, j -aligned segment e is cut by a spare if and only if it is cut by terminal i or j. Furthermore, 
for I G {i, j}, the following conditions are all necessary for segment e to be cut by terminal i: 

(1) pi G e t 

(2) For all terminals h preceding £, p^ > mine/,,. 

(3) Terminal £ is not last in the order 

For the following, let e be an i, j-aligned segment. For probability distributions P on spares, one can 
obtain bounds on Tk(P, e) by using Conditions 1-3 above. For example, we can restrict our attention to 
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Condition 1: If pi and pj are uniformly distributed over [0,1], Condition 1 holds for terminal i with prob- 
ability |ej| = |e|, and similarly for terminal j. Thus, by linearity of expectation, the expected number of 
times e is cut is at most 2\e\. 

Next, consider adding Condition 3. Suppose that the ordering of terminals is random, meaning that i is 
last with probability l/k. The probability that e is cut by i becomes (1 — l/fc)|e|, so Tfc(P, e) < (2 — 2/k). 
Thus, uniformly random pfS and a random ordering gives a performance guarantee of of 2 — 2/k, matching 
the bound of Dahlhous et al. [ 5 1. 

To improve these bounds, one must use Condition 2. Calinescu et al. choose a spare by selecting p 
uniformly at random in [0, 1], setting pi = p for each terminal I, and slicing off terminals in random order. 
Conditions 1 and 3 again derive a density bound of 1 — l/k for terminal i and j. Calinescu et al. improve 
this analysis as follows. Suppose that the edge e is farther from j than from i. We will argue that the density 
contribution from terminal j (i.e. the contribution to Tk(P, e) due to terminal j cutting e) is only 1/2. The 
point is that if p is such that j potentially cuts e, i.e. p G ej, and if i (which is closer to e) precedes j in the 
random slice ordering, then i will capture all of e and prevent j from cutting it. 

Formally, we argue as follows. Without loss of generality, assume that min ej > max ej (note that any 
i, j aligned edge can be split in two with one part closer to i and one part closer to j, and our assumption 
then applies to each part separately). As was argued before, the contribution of terminal i is at most I — l/k. 
On the other hand, with probability 1/2, i precedes j. If so, since pi = pj, Condition 1 — that pj € ej — 
contradicts Condition 2 for i — that pi > minej. Thus e can only be cut by terminal j if j precedes i, in 
which case by Condition 1, the density contribution from j is 1. Thus the density contribution from terminal 
j is 1/2, leading to a total density of 3/2 — l/k. 

To improve on the 3/2 bound, we made stronger use of Condition 2. The analysis of Calinescu et al. only 
considers that a segment may be captured by the two terminals with which it is aligned. We derive stronger 
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results by observing that other terminals may capture the segment as well. To do so, we had to change the cut 
distribution as well as the analysis. It can be shown that no distribution that holds all pi equal can do better 
in the limit than the 3/2 factor of Calinescu et al. . For independent, uniformly distributed pi we also get the 
3/2 factor. The 3/2 factor can be improved somewhat in the limit by using a non-uniform distribution on each 
(independent) pi. However, the best cutting schemes we have found are based on combining dependence 
and non-uniformity. One such scheme for 3-way cut gives us a bound of 12/11, which is optimal over all 
schemes for 3-way cut. Another scheme gives us a bound of 1.3438 that holds for any number of terminals. 
This latter scheme is designed for large k; optimizing it for smaller k gives better bounds. 

2.4 Additional Observations 

What is the best embedding? Perhaps the first natural question to ask is whether the embedding chosen 
by Calinescu et al. is the best possible. 

Lemma 2.6 Among all relaxations based on embeddings in the simplex that minimize some norm ( without 
adding other constraints) the L\ norm has the smallest possible integrality gap. 

Proof: We show that the L\ norm maximizes the measured volume of any embedded graph; thus it mini- 
mizes the integrality gap. 

Suppose we use some norm || • ||' instead of the the (scaled) L\ norm | • |. If the norm provides a 
relaxation, the distance between simplex vertices must be at most one — that is, for any edge e connecting 
simplex endpoints, we have Hell' < 1 = |e|. Consider some embedded edge e. As discussed in Section l2~2l 
we know that under the L\ norm it is connected by a path of aligned edges e\, . . . , e r such that |e| = J2 \ e i\- 
Since any norm-based distance measure is translation invariant and proportion preserving, this implies that 
for each a (which is a scaled, translated version of an edge of the simplex) we have ||ej||' < \e{\. It follows 
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from the triangle inequality that 



||e||' < X^II^H' 

< EN 

Since || • ||' assigns no greater a length to every embedded edge than | • |, it also assigns no greater a volume 
to any embedded graph. Thus, its integrality gap is no better than that induced by the L\ norm. □ 

Symmetry. A second observation is that there is no benefit in trying to identify a "good terminal order" in 
which to cut up the simplex. 

Lemma 2.7 There is a spare cutting scheme whose maximum density is minimum among all spare cutting 
schemes and that has the following form: 

1. choose slice distances (di, . . . , dk-i)from some probability distribution 

2. apply the slice distances (in order) to a uniform random permutation of the terminals 

An analogous "order independence" statement holds for the best possible (possibly non-spare) algorithm. 
Proof: For any cutting scheme P, let P' the the corresponding "symmetrized" cutting scheme, i.e. 

Pr[P' cuts corner 1 at distance d\, then cuts corner 2 at distance di, etc.] 
= £f J2a P r [-P cuts corner cr(l) at distance d\, then. . . ] 

where a varies over all permutations of 1, . . . , k. 
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For a line segment e and a permutation a, let <r(e) denote the line segment obtained by permuting by a 
the coordinates of the start point and end point of e. 
Then for any line segment e, 



density of P' on e = — ^ density of P on cr(e) 



Let r be the maximum density of P. Then for any a, the density of P on <r(e) is at most r. The density 
of P' on e is thus the convex combination of values all of which are at most r, so is in turn at most r. Thus 
the density of P' is no more than that of P. q 

The above lemma shows that there is no worst-case benefit to considering specific terminal ordering. 
The duality argument of Section 12.11 carries over to show that a spare with optimum expected maximum 
density can be specified simply as a distribution over slicing distances, without reference to an input graph 
embedding. 

3 Our Computational Study 

In this section we describe some computational experiments we carried out to help us understand the behav- 
ior of the geometric embedding. These experiments also yielded the best rounding schemes so far for the 4- 
and 5- terminal cut problems. One need not read this section in order to understand the following ones. 

As discussed above, our goal was to find a distribution over cuts of the /c-simplex that minimized the 
density of any segment in the simplex. This problem can be formulated as an infinite dimensional linear 
program, with one variable per cut of the simplex, corresponding to the probability that that cut is chosen, 
and one constraint for every (aligned, infmitesimaHy small) line segment inside the simplex, which measures 
the expected number of times the chosen cut will cut that segment. Of course, it is not tractable to solve the 
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infinite LP computationally, but we expected that discretized versions of it would be informative. 

We applied this approach in two distinct ways. For the 3-terminal case, we devised an LP that exploited 
the planarity of the 3-terminal relaxation to home in on a "worst case" embedded graph. By examining this 
graph, we were able to deduce requirements for the optimal algorithm, which led to its identification. For the 
general case, we devised an LP whose solutions are (provable) upper bounds on the performance of certain 
rounding algorithms. We solved this LP for small numbers of terminals (3-9), deriving algorithms with 
(computer aided) proofs of the best known performance ratios for these problems. The solution suggested 
certain properties that appear to hold in the "optimal" rounding scheme; we used these suggestions in our 
development of (analytic) solutions for arbitrary numbers of terminals. 

3.1 The three-terminal case 

For the 3-terminal problem we exploited planarity. The 3-simplex can be viewed as a triangle in the plane. 
We discretized the linear program by defining a triangular mesh over the simplex and considering only 
edges of the mesh instead of all line segments in the simplex. A (rather coarse) example mesh is shown in 
Figure [2 Note that we have augmented the triangular mesh with rays starting at the corners of the simplex 
and heading out to infinity. 

We used the planarity of the 3-simplex to simplify our LP formulation. The planar dual of the augmented 
mesh is shown in Figure |2] Note that because of the augmentation, the dual has three auxiliary nodes A, B, 
and C corresponding to infinite regions of the primal, one node for each side of the simplex. 

Any minimal 3-way cut of the mesh corresponds to a collection of two or three paths (representing the 
boundary of the cut) through the planar dual of the augmented mesh. Specifically, the cut corresponds to 
either 

1. two paths whose endpoints are the nodes A, B, and C (illustrated in Figure^, or 
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Figure 1 : A triangular mesh used in the linear program for the three-terminal case. The mesh is augmented 
with rays going from the corners of the simplex to infinity. 




Figure 2: The planar dual is shown in bold. The nodes A, B, and C correspond to the three infinite regions 
of the primal. 
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c 

Figure 3: An example of a 3- way cut corresponding to a pair of paths. 




2. three paths all originating at some interior dual node and terminating at the nodes A, B, and C 
(illustrated in Figure|4]i. 

Given an assignment of weights to the edges of the mesh, the weight of the minimum 3-way cut of the first 
type is 

min sum of distances in dual graph from X to each of the other auxiliary nodes 

Xe{A,B,C} 

and the minimum 3-way cut of the second type is 

min sum of distances in dual graph from x to A, B, and C 

x ordinary node of dual 

To find an embedded graph that is "worst-case" (up to the discretization), we solve a linear program 
that has a weight variable w uv for each edge uv of the triangular mesh (not including the rays). The linear 
program also has a distance variable d xy for every pair x, y of nodes of the dual (including auxiliary nodes). 
The objective is to minimize the total weight J2 e x e subject to the condition that every 3-way cut has value 
at least 1 . This condition can be expressed by a collection of constraints on distances through the dual graph. 



min w uv s.t. 



dxA + d x s + d x c < 1 for each ordinary dual node x 

dAB + d A c > 1 
dBA + d B c > 1 
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dcA + d C B > 1 

d X z < d xy + w yz for each ordinary dual edge yz 
d X x = for each dual node x 

w uv > for each edge uv 

d xy > for each dual edge xy 

Using the above linear program, we first deduced the general form of the dual solution, giving us the 
lower bound for k = 3. From this we deduced the necessary structure of any optimal primal solution (using 
complementary slackness conditions), including the important idea of "ball cuts" versus "corner cuts" which 
we will discuss in the following sections. 

3.2 The general case 

In the general case, the lack of a planar embedding prevented us from exploiting nice properties of its cuts; 
we were faced with the problem of enumerating cuts as well as edges. Based on the work of Calinescu et 
al. and our own results for the optimal 3-terminal solution, we decided to limit our exploration to spares as 
discussed above. 

There is still an infinite space of possible spares, so we discretized our problem. Fix an integer grid size 
N. A discrete spare is described by a vector (qi, . . . , qu-i) where each qi is an integer in the range [0, N—l]. 
Given such a vector, we choose a random spare by setting di uniformly in the range [qi/N, (q^ + 1)/N]. 
This defines a probability distribution on spares. We now define a linear program to search for a probability 
distribution over all discrete spares (which induces a probability distribution over all spares). We define a 
variable for each discrete spare, which reflects the probability of choosing that discrete spare, and provide 
constraints that upper bound the density of every possible segment under this probability distribution. We 

18 



then aim to minimize the largest of these densities. 

There are infinitely many segments, but we define a finite set of constraints that allow us to upper bound 
the density of all of them, as follows. The slices at distances q/N (q = 1,2, . . . , N — 1) for each terminal 
partition the simplex into cells 

{(x u ...,x k ) : qi /N < Xi < (zi + 1)/N} 

For a given distribution on the discrete spares, we can compute a (linear) upper bound on the density induced 
on any segment with a given alignment within a cell, and specify one constraint saying that this upper bound 
should be small. Since the cells are small, we expect all segments with a given alignment to have roughly 
the same density under our cutting scheme, so we hope that the upper bound is reasonably tight. With this 
simplification, the number of constraints is bounded by the number of cells times the number of segment 
alignments per cell, which is at most k 2 N k . 

We determine the upper bound for a cell as follows. For any discrete cut, the slices generated from it 
will fall into one of three categories. If the i th coordinate of the discrete cut is different from that of the cell, 
then the i th slice will not pass through that cell: depending on whether the coordinate is larger or smaller it 
will either capture the entire cell or none of the cell. If the i th coordinates are the same, then the slice might 
pass through the cell; we can use the fact that the slice is uniformly distributed over a range to determine its 
density contribution. 

An i, j-aligned segment can only be cut if the slices for terminal i or j go through its cell (and no earlier 
slice captures the entire cell). If only one of the two slices goes through the terminal then its contribution to 
a segment's density is at most N (the length of the segment divided by the width of the cell). If both slices 
go through the cell, their contribution is at most 2N. We ignore the fact that different slices within the cell 
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might capture the segment before it can be cut, thus introducing some slack in our upper bound. 

In view of Lemma |2~71 we need only represent cutting schemes that consider terminals in random or- 
der. Recall that each assignment of integers in [0, N — 1] to q%, . . . , q k -i defines a discrete spare. We let 
x qi,-,qk-i b e tne variable assigning a probability to this spare in the cutting scheme. For any permuta- 
tion a of {1, ... , k}, the probability of applying this spare to the sequence of terminals a(l), . . . , a(k) is 

fc! x <?iv,<Jfc-i- 

In order to present the linear program, we require a bit of notation. For qi, . . . , qk-i, ai, • • • , Qfe € 
{0, . . . , N — 1} (the qgs representing a discrete spare and the a/s representing a cell), for a permutation a 
of {1, ... , k}, and for i £ {1, . . . , k}, define 

{0 if 3m < cr~ 1 (i) : q m > a a(jri ) 
1 otherwise 

The value of ff is if some slice earlier than the slice for terminal i captures the cell defined by a\, . . . , a^. 
Define <5(m, n) to be 1 if m = n and zero otherwise. 
The linear program minimizes r subject to the following constraints. 

x qi, --Ak-i = 1 (1) 

qi,—,qk-l 



T|5Z (f?((qi,---,qk-i),( a ii---, a k))5(Qa-Hi),ai)N (2) 



qi,—,qk-i a 



+ fj({qi,---,qk-i),(ai,---,a, k ))S(q a -i ( j- ) ,a j )N)x qit ,,, t g k _ 1 <r (3) 



We can exploit symmetry to further reduce the number of constraints we consider. Since by assumption 
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our spare slices terminals in random order, two segments that are identical under permutation of coordinates 
will have the same densities, so we need consider only one of them. Thus, we restrict our constraints to 1, 2- 
aligned segments in which the remaining coordinates are in nondecreasing order. For each such segment and 
cell, we use one constraint to measure the average density induced by a given spare over all permutations of 
the terminal orders. 

3.3 LP Results 

Exploiting symmetry as discussed above, we were able to solve relatively fine discretizations of the prob- 
lem. We wrote a simple program to generate the linear programs automatically, and used CPLEX to solve 
them. While it is difficult to "prove" programs correct, our computations did converge to the correct 12/11 
approximation ratio for the 3-terminal case. 

We give our results below in tabular form. We derived improved bounds for 4-9 terminals. Note that 
these programs optimize a proven upper bound on the approximation ratio; thus,under the assumption that 
the programs were correct, these numbers are proven upper bounds. In fact, since the programs output a 
particular distribution over discrete cuts, their performance ratio could be proven analytically via a tedious 
case analysis on each cell of the discretize grid (which we have not performed). 



k 


grid 
size 


corner cut 
probability 


bound 


3/2 - 1/k 


3 


90 


.284 


1.0941 


1.167 


4 


36 


.289 


1.1539 


1.250 


5 


18 


.314 


1.2161 


1.300 


6 


12 


.376 


1.2714 


1.333 


7 


9 


.397 


1.3200 


1.357 


8 


6 


.414 


1.3322 


1.375 



Our experiments also revealed one interesting fact: in all cases, the optimum cut distribution made use 
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of "corner cuts." That is, the output distribution had the following form: with some probability, place each 
slice at a distance chosen uniformly between and 1/3 from its terminal; otherwise, use a (joint) distribution 
that places every slice at distance greater than 1/3 from its terminal. 

Adding constraints that forced the corner cuts to operate over a range other than 1/3 of the way from 
the terminals worsened the computed performance ratio, hinting that perhaps the optimal algorithm uses 
corners of size exactly 1/3. This result is consistent with the optimal 3-terminal algorithm, but inconsistent 
with the corner cut placement in the analytical solution for higher k that we give later. We may be observing 
a misleading artifact of working with a small discretized problem, or we may be missing something in our 
analytic solution. 

4 Upper Bound for k = 3 

Our analytic upper bound of 12/11 for k = 3 comes from a new cutting scheme that we call the ball/corner 
scheme. Though for simplicity we present a non-spare scheme, there is a similar scheme using spares that 
achieves the same bound. 

For k = 3, the simplex A can be viewed as a triangle in the plane, which simplifies our pictures. How- 
ever, we continue to use the original three-dimensional coordinate system to locate points in the simplex. 
Our cut of the simplex is determined by some lines and rays drawn through the triangle; we refer to them as 
boundaries. We will show that no segment has high density with respect to our random choice of boundaries. 

As illustrated in Figure|5j number the vertices of the simplex 1, 2 and 3. Let points a,b,...,f divide the 
edges in thirds, so that a-b-f-d-c-e-a is the hexagon in A with side length 2 1/3, with side c-d on the side 
of the triangle connecting terminals 2 and 3. Note that this hexagon is (a scaled version of) the unit ball for 
our distance metric. The points on the boundary of the hexagon are each at distance 1/3 from the hexagon's 
2 Remember that we measure length as half the L\ norm. 
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Figure 5 : This figure illustrates the cuts used for the case k = 3. The ball is contained within the dotted 
lines. The leftmost diagram shows how r might be chosen for the ball cut. The middle diagram shows one 
possible resulting ball cut (bold lines). The rightmost diagram shows a corner cut (bold lines). 

center.Outside the hex, we have a corner for each terminal i consisting of the points x with X{ > 2/3. 
4.1 The Ball/Corner Scheme 

The ball/corner scheme chooses a ball cut with probability 8/11, otherwise it chooses a corner cut. These 
two types of cuts are defined next. The scheme is illustrated in Figure |5] 

Ball cut: Choose a point r uniformly at random from either line a-c or line b-d. Consider the three lines 
A Xi=n (i = 1, 2, 3) parallel to the triangle's sides and passing through the point r. Each such line is divided 
at the point r into two rays. Thus we get six rays. Each side of the triangle intersects two of these rays. 
For each side, choose uniformly at random one of the two rays that hit it. This gives three rays; they form 
the boundary of the 3-way cut. (For a spare-based equivalent of this algorithm, we can choose (in random 
order) two of the three spares that pass through the chosen point r.) 

Corner Cut: Choose two terminals in {1, 2, 3}, and a value p € (2/3, 1], uniformly at random. For each 
of the two chosen terminals i, let U = A Xi=p . The two lines U form the boundaries of the 3-way cut. (Note 
that a corner cut is a spare.) 
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4.2 Analysis 

We first state two simple properties of the ball cut that we need to analyze the performance of the cutting 
scheme: 

Fact 4.1 Each of the 3 coordinates of the random point r is uniformly distributed in [0, 2/3]. 

Fact 4.2 Once r is chosen, each one of the six candidate rays connecting r to one side of the triangle is 
chosen with probability 1/2. 

Theorem 4.3 The maximum density of the ball/corner scheme is 12/11, so r| < 12/11. 

Proof: We show that the expected density of any segment e is at most |e| • 12/11. For the ball cuts, we use 
only the two facts claimed above. Since these two facts, as well as the corner cut scheme, are symmetric 
with respect to the three coordinates, it suffices to prove the claim only for a 1, 2-aligned segment e. Further, 
we may assume assume that e is entirely contained in either a corner or the hex; for otherwise, as discussed 
in Section im we can just split e into corresponding pieces, calculating the density for each piece separately. 
We will consider several cases, depending on where e is located. 

First, assume e is located entirely in the hex. Such a segment cannot be cut by a corner cut, so we need 
only consider the density when a ball cut is made and multiply by the probability of choosing a ball cut, 
namely 8/11. Assume a ball cut is made. Then e can only be cut by rays in A Xi=ri fori = 1,2. By Fact PTTl 
ri is uniformly distributed in [0, 2/3]. Hence, the probability that A Xi=rz goes through e is |e|/(2/3) since e 
is 1, 2-aligned. If A Xi=n touches e, it is at a single point. By Fact l4.2l the ray of A Xi=ri containing this point 
is picked for the cut with probability 1/2. Thus the expected number of times e is cut is yy • 2 • ^ • ^ = yy |e|. 

Exactly the same argument applies if the edge is in the corner closest to terminal 3. The ball cut con- 
tributes the same 12/11 density, while the corner cut contributes nothing (note that a 1, 2-aligned edge is 
parallel to the line A X3=r3 , so cannot be cut by it). 
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Finally, suppose segment e is in the corner closest to terminal 1 (a symmetric argument applies if e is 
in the corner closest to terminal 2). In this case, if a ball cut is made, the above analysis applies except 
that only the line A X2=r2 can cut e (the line A xi=ri never enters the corner), so the density contribution of 
the ball cut is halved to |e|yr. But the edge can also be cut by a corner cut. A corner cut is chosen with 
probability 3/11. When it is, two of the three terminals are chosen, so terminal 1 is chosen with probability 
2/3. If terminal 1 is chosen, then, since the cutting line near terminal 1 is of the form A xi= ±_ p , where p is 
chosen uniformly in [0, 1/3], the probability that the line cuts e is |e|/(l/3). Thus, the expected number of 
times that the edge e is cut (by a ball cut or corner cut) is |e|^j- + -jj ■ | • = [e| 

□ 



5 Lower Bound for k = 3 

Theorem 5.1 For k = 3, the minimum maximum density r| > 12/11. Hence, the integrality gap for the 
geometric relaxation is 12/11. 

Note that this theorem applies to all cutting schemes, not just spares. Thus, the scheme of the previous 
section is optimal. 

Proof: Fix N to be any positive integer. We construct an embedded weighted graph Gn with no 3-way 
cut of cost less than 12iV, but with an embedding of cost lliV + 1. This immediately demonstrates an 
integrality gap of Y1N/{\1N + 1). Furthermore, it implies that no cutting scheme has maximum density 
less than \2N/{1\N + 1), because by Lemma |2~T1 such a cutting scheme applied to Gn would yield a 
3-way cut with expected cost less than \2N , a contradiction. Since N is arbitrary, the result follows. Our 
construction (for N = 7) is shown in Figure |5| 
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Figure 6: The lower bound for k = 3 (here TV = 7). The paths from 2 to 3 are on the left. The entire graph 
is on the right. On the border, overlapping paths are drawn side-by-side for clarity, so line width represents 
edge cost. 

For any pair of distinct terminals i, j and number d € [0, 1], define embedded path p(i,j, d) as follows. 
Let i be the terminal in {1, 2, 3} — {i, j}, let a be the point on segment it at distance d from i, and let b be 
the point on segment jl at distance d from j. Then p(i, j, d) is the union of the three segments ia, ab, and 
bj. 

We form the graph from 9iV paths p(i,j, d) for < d < 2/3, where d is an integer multiple of l/(3iV). 
Although we describe the graph as a set of paths, technically it is a planar graph consisting of nodes and 
edges as follows: for every point in A whose coordinates are integer multiples of 1/(3 AT), there is a node in 
the graph embedded at that point; for every pair of nodes embedded 1 / (3N) units apart, G has an edge with 
cost equal to the number of paths that pass through both nodes. 

With this understanding, we now specify the graph. For each of the 3 distinct pairs of terminals 
there are 3iV paths. Of these paths, N run directly between the terminals; that is, there are N copies of 
p(i, j, 0). The remaining 2N paths are the paths p(i, j, mj (3N)) where m = 1,2,..., 2N. 

The total cost of the embedding is the total length of the paths. Since a path p(i, j, mj (3iV)) has length 
1 + m/(3N), a direct calculation shows that the total length of the paths is 3[iV + 1 + m/(3N)\ = 

1L/V + 1. 

Next we lower bound the cost of any 3-way cut. Since the graph is planar, any minimal 3-way cut 
corresponds either to a disconnected cut (meaning that the cut is the union of two disjoint 2-way cuts, each 
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separating some terminal from both other terminals), like our upper bound's corner cut, or a connected cut 
(meaning that the cut edges give, in the planar dual, three paths connected at some central node and going 
to the three sides of the triangle), like our upper bound's ball cut. 

Any 3-way cut must cut all of the 9N paths at least once. To finish the proof, we will argue that for 
either type of 3-way cut (connected or not), at least 3N paths are cut twice, so that the edges cut by the 
3-way cut cost at least 12N. This is easy to verify for a disconnected cut: a disconnected cut is the union of 
two 2-way cuts, so the 3N paths running between the two terminals that are cut off must be cut twice. 

Now consider any connected cut. In the planar dual of Gat, the connected cut corresponds to a central 
node and three paths from the node to each side of the triangle. Let x = {x\, 22, ^3) be any point inside 
the face of Gn corresponding to the central node. Consider a path p(i,j,d) such that d > Xf, where 
I 7^ i, j. That is, x is inside the cycle formed by the union of p(i, j, d) and p(i, j, 0). Then the path p(i,j, d) 
is cut at least twice by the connected cut. For fixed i and j, the number of such paths (with d > xi) is 
|_2iV— xe3N+l\ > 2N—3xeN. Thus, the total number of such paths is more than 6iV— 3(xi+X2+xs)N = 
3N. n 



6 Improvement for general k 

We now present an algorithm for an arbitrary number of terminals. While this algorithm seems unlikely 
to be the best possible, it improves on the previous best bound. As discussed in Section 12.31 the essential 
observation in this analysis is that many slices can capture an edge before it has a chance to be cut. 

Theorem 6.1 For all k, t£ < 1.3438. Moreover, there is a k-way cut randomized approximation algorithm 
with an approximation guarantee 0/1.3438. 
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Our bound improves on the Calinescu et al. bound of 1.5 — 2/k for all k > 14. For k < 14, we show that 
rt < 1.5 — 2/k by specializing the analysis for small k (see Subsection 16. II) . 

To prove the theorem, we will use a (spare) cutting scheme, that is, we choose k slicing thresholds pi, 
and apply the slices A Xi=Pi to a random permutation a of the terminals. We are going to apply either an 
independent cut (ICUT) or a corner cut: 

ICUT: each pi is chosen independently and uniformly in [0, 6/11]. 

Corner cut: all p^ are chosen equal to a single random p picked uniformly in [6/11,1]. 

We will apply ICUT with probability a = 0.667186 and apply a corner cut with the remaining probability. 

Before proving that the above spare achieves a maximum density below 1.3438, as in Theorem l6.ll let's 
first draw parallels to our scheme for k = 3. The corner cut is completely analogous, except that p is now 
chosen in the interval [6/11, 1] instead of [2/3, 1]. However, ICUT is very different from the ball cut. For 
example, the pi are now independent whereas they were highly dependent in the ball cut. The reader may 
wonder why we did not just generalize the ball/corner cut scheme. However, corner cuts are only meaningful 
for p = (1/2, 1] since this is the region in which the A Xi=p are disjoint. On the other hand, the maximal 
ball of the fc-simplex has center (1/k, . . . , 1/k) and radius 1/k. Already for k = 4, the simplex cannot 
be covered by a ball and corners, and for k large, the measure of the ball is vanishingly small. Hence, the 
concept of ball cuts is not really relevant for large k. 

To bound the cutting density of our scheme, we will bound the density of every segment. As justified in 
Section |2~21 we consider a segment of length e > 0, and let e approach zero. As in the ball/corner scheme, 
by symmetry we can assume without loss of generality that the segment is 1 , 2-aligned. 

Define df.(xi, . . . , x^) to be the density with which ICUT cuts a 1, 2-aligned segment of infinitesimal 
length located at x\, X2, • • • , x^. We will show: 
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Lemma 6.2 

( 2.014096 ifx u x 2 < 6/11 
d k (xi, ...,x k )< < 

I. 11/12 otherwise. 

By Lemma l6~2l this combined scheme gives a density of 2.014096a for non-corner segments (since 
they are cut only if ICUT is used, and then only with probability 2.014096), and a density of (ll/12)a + 
(11/5)(1 — a) for corner segments (combining their probabilities of being cut by the two schemes), for a 
maximum density of max{(2.014096)o, (11/12) a + (11/5)(1 - a)} < 1.3438, proving Theorem O To 
finish the proof of Theorem 16. II it remains only to prove Lemma |6~2"1 

ICUT's cumulative probability distribution function for any pi is F(z) = mm{(ll/6)z, 1}. The corre- 
sponding probability density function is 

{11/6 if* €[0,6/11] 
otherwise. 

Consider a 1, 2-aligned segment of length e with one endpoint fixed at x\, x 2 , ■ ■ ■ , Xk- As e goes to zero, the 
density of this segment goes to 




where the sum is over all k\ orderings of the terminals. This formula follows from Fact 12.51 The first term 
measures the probability that the segment is cut by terminal 1 , which happens if the slice for terminal 1 goes 
through the segment while all slices preceding terminal 1 in the ordering fail to capture the segment. The 
second term similarly measures the probability that the segment is cut by terminal 2. Considering the slices 
for terminals other than 1 and 2 is the crucial element in improving the density bound of 3/2 for large k. The 
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formula assumes that F is continuous around each Xi and that F' is continuous in an open region around 
x\ and x 2 . The latter is not the case around 6/11. However, as discussed Section |2~71 we may assume 
that all segments e considered have been subdivided so that for each i = 1,2, either maxxj < 6/11 or 
minxj > 6/11. 

Note that d k {x\, ■ ■ ■ , Xi, 0, . . . , 0) = di(x\, . . . , xi) (provided i > 2), because Xj = implies terminal 
j cannot save the edge. Note also that d k is symmetric with respect to the variables Xj for i > 2. Define 

D k (xi,x 2 ) = max d k (x 1 ,x 2 , . . . ,x k ) 

x 3 ,...,x k 

C k (xi,x 2 ) = d k (x 1 , x 2 , c, ... ,c) 

where c = (1 — x\ — x 2 ) /{k — 2), 
-DooOi,£2) = lim D k (xi,x 2 ), 
C 0O (xi,x 2 ) = lim C k (xi,x 2 ). 

k—*oo 

In these definitions, (x%, x 2 , . . . , x k ) is required to lie in the fc-simplex. 

D k (xi,x 2 ) is the maximum density of any 1, 2-aligned infinitesimal segment with an endpoint whose 
first two coordinates are x\,x 2 . Note that the maximum is well-defined and achieved by some x%, ■ ■ ■ ,x k 
because the simplex is closed under limits. 

To understand ICUT, our first goal is to characterize D k . We consider C k as it is one candidate for D k . 

Lemma 6.3 D k (x\,x 2 ) < D k+ i(xi, x 2 ) for all k. 

Proof: max(4(xi, ... ,x k ) = maxd k+ x(xi, ... ,x k ,0) < maxd k+ i(x 1 , . . . ,x k ,x k+ i). □ 

Thus for fixed x\,x 2 , {D 2 (xi, x 2 ), D%(xi, x 2 ), . . .} is a nondecreasing sequence bounded from above 
(by 2). This implies that is well-defined. We will see later that is also well-defined. 
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Next we show that for fixed x\ and X2, the maximum for d k occurs at either the "central point" 

xi, X2, c, c, . . . , c or the "three-terminal" point xi, x 2 , 1 — x\ — X2, 0, . . . , 0. 
Lemma 6.4 

(C k (x 1 ,x 2 ) ifVi > 2 : Xi < 6/11 
d k (x!, . . . ,x k ) < < 

lc 3 (x 1 ,x 2 ) i/3i > 2 : Xi > 6/11. 

Proof: Fix x\ and X2- Let c = (1 — x\ — x 2 )/{k — 2). 

Claim 1: Among all X3, . . . , x k such that < X{ < 6/11 for all i > 2 (orai xi, X2, . . . , x k is in the 
simplex), the unique maximizer of d k (x±, X2, X3, . . . ,x k ) satisfies X3 = X4 = • • • = x k , so is equal to 
C k (xi, X2). Suppose for contradiction that some other such x 3 , x 4 , ... ,x k maximizes d k . Then x« < Xj for 
some i,j > 2. Considered just as a function of x« and Xj (holding the other coordinates fixed) 

d k (x u ...,x k ) =p + q[l-F{xi)]+r[l-F{xj)\ + s[l - F{ Xi )\ [1 - F{x 3 )} (5) 

where p, q, r and s are nonnegative and independent of Xj and X j. Furthermore q = r because d k is 
symmetric in Xj and Xj. Consider increasing x, and decreasing xj at equal rates. This maintains < 
Xj, Xj < 6/11 but increases d k at a rate proportional to 

q[F'( Xj ) - F'( Xi )] + s(F'( Xj )[l - F( Xi )] - F'( Xi )[l - F( Xj )]). 

This is positive because F'(z) = 11/6 for z < 6/11 and F(xj) > F(xi) (recall that x« < Xj < 6/11). This 
contradicts the choice of X3, . . . , x k . 

Claim 2: Among all X3, . . . , x k such that x« > 6/11 for some i > 2 (and Xi, . . . , x k is in the simplex), 
the unique maximizer of d k (xi,x 2 , x 3 , . . . , x k ) satisfies X{ = 1 — x\ — X2 and Xj = for j ^ i, so is 
equal to Cs(xi,x 2 ). Suppose for contradiction that some other such x 3 ,X4, . . . ,x k maximizes d k . Fix 
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some j > 2 such that < Xj < 6/11 < X{. Since by assumption Xj(?e6/ll, we have F{xi) = 1 and thus 
the expression ® reduces to p + r(l — F{xj)). If we increase X{ and decrease Xj at the same rate, the rate 
of increase in dk is rF'(xj) > 0, contradicting the choice of x%, . . . , 24. 

The two claims together prove the lemma. q 

Lemma 6.5 Fork > 4, Ck(xi,X2) < Ck+i(x\, x 2 ). 
Proof: 

C k (x 1 ,x 2 ) = d k (x 1 ,x 2 ,c, ... ,c) 

= 4+i(xi,x 2 ,c, . . . ,c,0) 
< C fe+ i(xi,x 2 ). 

Here c = (1 — xi — x 2 )/{k — 2). The last inequality follows from Lemma l6T4l (using c < 1/2 < 6/11). n 

An immediate corollary is that C 00 (xi, x 2 ) is well-defined and Ck{x\,x 2 ) < C 00 (x\,X2) for all k. 
Using this and Lemma|OJ to bound it suffices to bound C3 and C^. We begin with C^. 

Lemma 6.6 

r 2.014096 i/xi,x 2 < 6/11 
C 00 (xi,x 2 ) < < 

[ 11/12 otherwise. 

Proof: Fix x\ and x 2 . Our first goal is to derive a closed-form expression for Ck(x\, x 2 ) for any k. Fix 
for now and let X{ = c = (1 — x\ — x 2 ) j(k — 2) for i > 2. 

For j = 1, 2, let Sj denote the probability that the segment at (x\, x 2 , . . . x^) is not captured by a 
terminal other than j before terminal j's cut is made: 

n (i-n*))- 

° i:a{i)<a{j) 
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Then C k { Xl ,x 2 ) = SiF'(xi) + S 2 F'{x 2 ). 

We will derive a closed-form expression for Si (and by symmetry for S 2 ). Recall that x\ = c for i > 2. 
We thus rewrite 

Si = I £ [1 - ^(c)]- 1 [1 - F(* 2 )] + (l - ^) [1 - F(c)]«) . 

<j=0 

Here we condition on q, the number of z such that a(i) < a(l). Note that q is uniform in {0, 1, . . . , k — 1} 
while is the probability that a (2) < cr(l), given g. A change of variables and rewriting give 

* = ( 1+i i^r 1 ) S 11 ^ - ^) S 2 ^- 

<7=0 <j=0 

Now we let k — > oo. The two sums above have standard closed forms that tend respectively to 

[l-e- a ]a- 1 and [1 - (1 + a)e~ a ] a~ 2 , 
where a = linifc^oo k F(c) = (1 — x\ — x 2 )F'(0). Thus, as k — > oo, 

Si -f [1 - e^cT 1 - F(x 2 )[l - (1 + a)e~ a ] a" 2 . 

Of course S 2 tends to the above with xi replacing x 2 . This gives us our closed-form expression for 

C 00 (xi,x 2 ): 

Coo(xi,x 2 ) = [F\xi) + F\x 2 )\ x - [F'{xi)F(x 2 )+F'(x 2 )F{xi)] x 1 ~( 1 + °) e (6) 

where a = (1 — xi — x 2 )F'(0). 
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The above equality holds for any suitably well-behaved F. Using this closed form and our particular 
choice of F, we now show the two desired bounds on C^. 



Case 1: x\, x 2 < 6/11. In this case a = 11/6(1 — x\ — x 2 ), F'(x{) = F'(x 2 ) = 11/6, and F{x{) + 
F(x 2 ) = ll/6(xi + x 2 ) = 11/6 - a. So © gives 



where a = 11/6(1 — x\ — x 2 ) so a € [0, 11/6]. Let C(a) = Coo(xi, x 2 ). In the rest of this case (Case 
1), we will prove that C(a) < 2.014096 for a G (0, 11/6). The cases a = and a = 11/6 follow by the 
continuity of C. The claim is "obvious" from a plot but the somewhat technical proof appears below. 

We show that C(a) is strictly concave for a G (0, 11/6). It therefore has a unique maximum at some 
a , where C'(a ) = 0. By substitution, C"(.294) > 0.00045 > and C"(.295) < -0.00009 < 0, so 
a G (.294, .295). Hence 




max C(a) < CY.295) - 0.001 • C"(.295) < 2.014096 

ae[0,H/6] 



To show C(a) is strictly concave, we show that C 



(a) is strictly negative. Now, 




11 7e~ a a 2 -18a- 4e~ a a 



11 6e- a a 3 + 22 - 22e 



—a 



a: 



3 




a: 



3 



and 
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(7 e~ a a 3 + 3 e~ a a 2 - 36 a - 30 e~ a a + 6 e~ a a 4 + 66 - 66 e~ a ). 



36a 4 
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To show that C"{a) is negative, it suffices to prove that 



D(a) = -7e~ a a 3 - 3 e~ a a 2 + 36 a + 30 e~ a a - 6 e" a a 4 - 66 + 66 e~ a 

is negative. By substitution, D(0) = and D(ll/6) = 0, so it suffices to show that D' has only one zero 
oi, D'(a) < for a < a\ and D'(a) > for a > a\. Here 

D'(a) = -17e- a a 3 -18e- a a 2 -36e- a a + 36 - 36e- a + 6e" a a 4 

and D"(o) = e" a a 2 (-6a 2 + 41a - 33). For a G (0, 11/6], L>" has only one zero a 2 = 41 ~ 1 ^ w 0.93 
and D"(a) < for a < a 2 and D"{a) > for a > a2. That is, D' is first decreasing and then increasing. 
Since D'(0) = and D'(ll/6) > 4.108 > it follows that D' has only one zero ai for a G (0, 11/6]. 

Case 2: xiorx2 > 6/11. Assume > 6/11 (the case x 2 > 6/11 is symmetric). In this case, F'(x\) = 
and F(x\) = 1, so we get 



„ . , 11 l-e~ a 11 l-(l + a)e- a 

C 00 (xi,x 2 ) = t. 

6 a 6 a z 



As before, let (7(a) = C 00 (x 1 ,x 2 ). We will prove that (7(a) < 11/12 for a G [0,11/6]. First, lim a ^ C{a) 
11/12, so (7(a) < 11/12 follows if we can show that C'{a) < for a G (0, 11/6]. We have 



C'(a) = ii(-a - e~ a a + 2-2 e - "). 
6a J 



Define E(a) = -a - e~ a a + 2-2 e" a . Since ^ > for a > 0, (7(a) < if and only if E(a) < 0. Since 
E(0) = 0, we can infer E(a) < if E'(a) < for all a G (0, 11/6]. We have E'{a) = -1 + e - a (a + 1). 
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Note that E'(0) = 0, so E'(a) < follows if E"(a) < for a € (0, 11/6]. We have E"(a) = -e~ a a, so 
E"(a) < 0. We conclude that Coo (xi,x 2 ) < 11/6 if xi > 6/11. □ 

Lemmas l6.4l thiough l6.6| prove that, for x such that Xi < 6/11 for all « > 2, 

d k (x x , . . . ,x k ) < C 00 (x 1 ,x 2 ) 

( 2.012096 ifxx,x 2 < 6/11 

< 

[ 11/12 otherwise. 
The remaining case is when xi > 6/11 for some i > 2. In this case by Lemma l6~4l 

d k (xi, ...,x k )< C 3 (xi,x 2 ) = d 3 (xi,x 2 , 1 - xi - x 2 ) 
and x\ + x 2 < 5/11. Thus, to finish the proof of the theorem, it suffices to show the following lemma. 



Lemma 6.7 Ifx 1 + x 2 < 5/11, then C 3 (xi,x 2 ) < 11/6 < 2.012096. 

Proof: Let x 3 = 1 — x\ — x 2 > 6/11. 

ThenF(aj 3 ) = 1 while F( Xl ) = ll/6xi,F(x 2 ) = ll/6x 2 ,and F'( Xl ) = F'(x 2 ) = 11/6. 

By inspection of &,C 3 (x u x 2 ) = d 3 (x u x 2 ,x 3 ) = (1/6) (11/6) (6 - 11/6 (x x + x 2 )) < 11/6. Q 

This proves Lemma l6~2l 
6.1 Improvements for small values of k 

For particular values of k it is possible to refine the analysis in the proof of Theorem |6. II to get improved 
bounds. In this case it is useful to modify the algorithm so that it only uses k — 1 cuts instead of k. In 
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particular, we do not use the cut for the terminal j with a(j) = k. The analysis for this modified algorithm 
goes similarly, with our definitions appropriately modified to reflect that we are using k — 1 instead of k 
cuts. 

Then, instead of passing to the limit, Ck(x\, x^) can be evaluated directly. Following this approach we 
obtained the following performance guarantees for particular k: 



K 


corner 
placement 


ICUT 
probability 


bound 


J 


.041 


.O/J 


1.131 


A 

4 


.oU/ 


.OOJ 


i.i sy 


c 

J 


. JOO 




1 OTX 


6 


.576 


.659 


1.244 


7 


.565 


.657 


1.258 


8 


.557 


.656 


1.269 


9 


.557 


.659 


1.277 


10 


.557 


.661 


1.284 


12 


.554 


.661 


1.293 


20 


.554 


.666 


1.314 


35 


.550 


.666 


1.327 



"Corner placement" is the placement of the corner (analogous to 6/11) and "ICUT probability" is the prob- 
ability of choosing ICUT. These parameters were chosen to try to minimize the resulting bound on the 
performance ratio, shown under "bound". These numbers are approximate; the ratios were evaluated nu- 
merically without formal verification. 
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7 Conclusion 



We have provided a better analysis of an embedding relaxation for multiway cut. We have exactly deter- 
mined the integrality gap for the 3-terminal problem, and given an approximation algorithm achieving the 
bound. For larger values of k, we have defined a class of cutting schemes called spares that, through a com- 
bination of nonuniform and dependent rounding, provide better approximation ratios than the best previous 
schemes. However, the question of the exact integrality gap remains open. 
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8 Appendix 

Proof of Theorem \2.3\ 

We show that there necessarily exists a cutting scheme whose maximum density equals the integrality 
gap of the relaxation. The basis of this proof is that the two quantities are solutions to dual linear programs. 
Although the linear programs are infinite, we show they have no duality gap. 

Interestingly, most of the proof holds in the following more general setting: We have a non-negative 
real vector x representing a relaxed solution to some problem. There is a set S of allowable solutions (also 
nonnegative real vectors), and we want to round x to some solution y € S. The method for rounding x is 
represented by a randomized rounding scheme, which is simply a probability distribution P on S. 

We assume that the cost of x is given by w ■ x = J2i WiXi for some nonnegative weight function w, and 
likewise the cost of any y € S is w • y. We want to choose a single rounding scheme P that has a good 
performance ratio against all possible weight functions w. (We will see that this is analogous to choosing a 
single rounding scheme of the simplex that has good performance ratio against all embedded graphs.) 

Define the performance ratio of a rounding scheme P to be 

E yEP [w ■ y] 

sup — . 

W w ■ X 

(This corresponds to the maximum density of a cutting scheme.) Here the notation E y€ p[w ■ y] signifies the 
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expectation over y € S chosen according to the probability distribution P. 



Define the integrality gap to be 



inf. 



y&S w-y 



sup 



w 



w ■ X 



the worst-case ratio of the minimum-cost of any true solution to the cost of the relaxed solution. 



Lemma 8.1 The performance ratio can be reformulated as 



sup 



E yeP [w ■ y] 



= sup 



w 



w ■ X 



In the case of k-cut, each index i corresponds to a "seglet" (edge of the divided simplex), and w cor- 
responds to an embedded graph (wi corresponds to the number, or total weight, of edges embedded along 
seglet i). In that case the lemma says that, to check the performance guarantee of a rounding scheme, it 
suffices to check it for each seglet. 

The proof is similar to the observation regarding two-player zero-sum matrix games that once one player 
has fixed their mixed strategy, the other player has an optimal mixed strategy that is pure. 

Proof: Clearly the left-hand side is greater than or equal to the right-hand side (take w to be any of the % 
unit vectors with a single coordinate equal to 1 and the others 0). To finish we show that the left-hand side 
is at most the right hand side for any fixed w. 

Fix w. Let A equal the right-hand side above. Note that by linearity of expectation, the left-hand side is 



w ■ X 
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By the definition of A, E ye p[yi] < Xxi, so the quantity above is at most 



W ■ X 



A. 



□ 



Theorem 8.2 If S and the dimension of x are finite, then there exists a rounding scheme P whose perfor- 
mance ratio equals the integrality gap. 

Proof: Let P y be the probability that we choose solution y. Choosing an optimum rounding scheme is 
equivalent to the following linear program: 



mimmizep r 



subject to < 



EyesPyf, < r (V.) 



Eugs 1 p y 



P y > (VyeS). 



The dual of this program is 



maximize^ A 



subject to < 



Ei ft = 1 



qi > (Vi). 
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By the change of variables Wi = qi/xi, this is equivalent to 



maximize^ A 



Ei W iVi 



> A 



(Vy € S) 



subject to < 



Ej WiXi 



= 1 



> 



(Vy G 5). 



But it is easy to verify that this is equivalent to the problem of choosing a weight function w to achieve the 
integrality gap. If x and the vectors in S are finite-dimensional and S is finite, then strong duality implies 
that the linear programs have equal values. rj 

Next we describe how this relates to the rest of the paper. Define edge set £ to be all edges (pairs of 
points) in the simplex. Define x e = |e| for e G £. Define S to contain the characteristic vectors y(C) of 
fe-way cuts C of the simplex: y e (C) = 1 if e is cut by C and otherwise. An embedding of a particular 
weighted graph G in the simplex then corresponds to a particular weight function w(G), where w e (G) 
equals the total weight of edges embedded on simplex edge e. With this x, S, and interpretation of w, we 
have the following correspondences: 
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general setting 


fc-way cut 


w 


<-> embedded G 


w • X 


= vol(G) 


yeS 


<-> A;-way cut of G 


w ■ y 


= cost of fc-way cut C of G 


rounding scheme P 


<-> cutting scheme P 


inf p performance ratio of P 


— T k 


integrality gap 


= integrality gap 



Thus, if S and the dimension of x were finite, we could conclude by the theorem that the t% equals the 
integrality gap of the relaxation. This yields an immediate corollary: 

Corollary 8.3 For input graph instances of any bounded size, there is a rounding scheme whose perfor- 
mance is equal to the integrality gap of the relaxation. 

Proof: Consider the set of all graphs whose size is bounded by some quantity. Each has a bounded number 
of vertices, and the (rational) weights on the graph are also of bounded size. The linear programming relax- 
ation thus also has bounded size. It follows that any vertex solution to the linear programming relaxation, 
which assigns embedding coordinates to all the vertices, has bounded size — meaning that the coordinates 
are rational numbers of bounded size. The set of coordinates at which vertices might be located in an optimal 
embedding therefore forms a discrete grid within the simplex, with a finite number of points. The embedded 
edges connect these points, so there is a finite number of embedded edges in the optimal solution. An input 
instance is determined entirely by the weights assigned to these edges, so has finite dimension. Similarly, for 
the purposes of rounding we need only consider fc-way partitions of the finitely many grid vertices. There 
are only finitely many of these partitions. 
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Since the dimension of inputs and the number of output solutions is finite, the previous theorem applies 
and shows that there is a rounding scheme with performance equal to the integrality gap of the embedding. 

□ 

If we want a rounding scheme that works (uniformly) for graphs of arbitrary size, we have to work 
somewhat harder. Wee show that the desired result follows as a limiting case of Theorem l8.2l 

In what follows, we restrict the A;-way cut problem to various particular subsets E of the edge set £ of the 
simplex. A fe-way cut of E is defined to be a subset C of E such that E — C contains no terminal-to-terminal 
path By the minimal maximum density with respect to E we mean 

t(E) = inf sup Pr[P cuts e]/|e|, 

p e&E 

where P ranges over cutting schemes of E (probability distributions over fc-way cuts of E). By the integral- 
ity gap with respect to E we mean 

gap(E)= S upinf ^ eeC 7 ( . e , ) | , 

where w ranges over weight functions on E with finite support. 

Our goal is to show r(£) = gap(£). Since we know gap(£) < r(£), it suffices to show r{£) < gap(f ). 

Define £ to contain those edges in £ with rational endpoints. (In fact, any countable dense set will do 
in place of the rationals.) Fix any enumeration of the edges in £. Let set £ n contain the first n edges in the 
enumeration. It suffices to show the following equality: 

swpT(£ n ) = t(£), (7) 
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because by Theorem l8.2l we know r(£ n ) = gap(£ n ), and clearly gap(£ n ) < gap(£), so combining with Q 
proves the theorem via 

t{£) = supr(£ n ) = supgap(£„) < gap(f ). 

n n 

Thus, to prove Theorem l2.3l we need only prove ©. In the remainder of the section we prove it as follows: 
We first show we can extend any sequence of cutting schemes, one for each £ n , to a single good cutting 
scheme for their union £ (edges with rational endpoints). This shows sup n r(£ n ) = r(f). We then show 
that we can extend the cutting scheme on £ to a cutting scheme on £ (all edges). This shows r(£) = t(£). 

One measure-theoretic issue that we must first address is how we formally define a probability distribu- 
tion P on our infinite sets £ and £ . For this we use Kolmogorov's Existence Theorem 13 ch. 7]. 

Before we explain, we introduce some terminology. Note that any cutting scheme P of a set F C £ 
induces a cutting scheme Pg on each subset E C F by restriction. We say that any P from which Pe can 
be so obtained is consistent with Pe- 

Kolmogorov's Existence Theorem implies the following. Fix any F C £. Consider a family of cutting 
schemes (Pe : E C F, E finite) of the finite subsets of £. If this family is consistent, meaning that whenever 
two cutting schemes Pe and Pe> in the family satisfy E C E', the second scheme is consistent with the first 
one, then there exists a single cutting scheme P of F, with P consistent with each Pe in the family. 

In this section, to describe any cutting scheme on £ (or £), we specify the consistent family of cutting 
schemes it induces on the finite subsets. The following lemma gives a useful condition for the existence of 
such a family. We will use it twice. 

Lemma 8.4 Let F be any subset of £. Suppose there exists a sequence of cutting schemes (Q^\ Q^ 2 \ ■ ■ •} 
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that converges in the following sense: for each finite E C F and each k-way cut C of E, the limit 

lim QP(C) (8) 

is well defined. Define Pe to be the cutting scheme of E that, for each k-way cut C of E, chooses C with 
probability in 0, i.e. P E {C) := lim 

n^oo Qe ^ {C)- Then (Pe '■ E C F, E finite) is a consistent family 
of cutting schemes, so by Kolmogorov's Existence Theorem there exists a cutting scheme P of F consistent 
with each Pe- 

Proof: We need to verify the following: 

1. Each Pe is a cutting scheme: J2c Pe{C) = 1 with each Pe{C) > 0, where C ranges over all cuts of 
E. 

2. For Pe and Pe> with E C E' both finite, Pe 1 is consistent with Pe- That is, for any cut C of E, 



p e (C) = J2Pe>(C) 

a 



where C ranges over all /c-way cuts of E' such that C n E = C. 

In each case, the desired property holds for the cutting scheme induced on E (and/or E') by for large 
enough m. For example, to verify the second property, use 

P E (C)-J2 p E'(C) = limQ^(C) - ^UxnQ^C') 

a n c 



lim 

n 

0. 



a 
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The first equality is by definition of Pg, the second is because the finite sum of the limits is the limit of the 
sums, and the last is because, for n large enough, the term inside the limit is well defined and necessarily 
(simply because is a cutting scheme, so the induced cutting schemes and Q^) are necessarily 
consistent). The first property follows similarly. q 

We need one last "utility" lemma. It will help us construct a sequence of cutting schemes (QW) that 
converges in the sense needed for Lemma lOl 

Lemma 8.5 Consider a countable collection C of countable sequences of real numbers, where each se- 
quence q = (gW , g( 2 ) , . . .) in C lies in some finite interval I s . 

Then there exists an infinite index set IC {1,2,...} such that for each sequence q € C, the limit 

lim q^ 

is well defined. 

Proof: The proof is a "dovetailing" variation of the standard proof that any sequence in a compact set 
contains an infinite convergent subsequence. 

Order the sequences arbitrarily and consider the infinite mauix whose jth row is the jth sequence. 
Associate with each sequence q an interval I' q , initially I q . Visit the sequences in the standard dovetailing 
order; that is, visit the jth sequence for each j = 1, 1, 2, 1, 2, 3, 1, 2, 3, 4, . . ., so that each sequence is visited 
infinitely often. 

While visiting a sequence q, narrow the associated interval to either its upper or lower half, and then 
delete from the matrix all the columns i such that q^ is no longer in the associated interval. Further, make 
the choice of upper or lower half so that infinitely many columns remain undeleted. (This is possible because 



47 



each of the infinitely many gW's lying in a column that was previously not deleted lies in one of the two 
halves.) 

To complete the construction, define i{n) to be the smallest index larger than n of any column that 
remains undeleted after the nth step of the construction, and take X = {i(n) : n = 1, 2, . . .}. It is easy to 
verify that for each q G C, the subsequence (gW : j gl) converges. q 

Now we can begin our two-step proof that t{£) = t{£) = sup„ r(£ n ). First we show that if each £ n 
has a good cutting scheme, then there is an equally good cutting scheme of £. 

Lemma 8.6 There exists a cutting scheme P of £ with maximum density sup n r(£ n ). Thus, t(£) = 
sup n r(£ n ) 

Proof: Let r sup = sup ra r{£ n ). Say a sequence of cutting schemes (P^) is good if each P^ is a cutting 
scheme of £ n and has maximum density at most r sup with respect to £ n . By the definition of r sup , there 
exists a good sequence of cutting schemes (Q^ : n = 1, 2, . . .}. 

Consider each pair (E, C) where E is a finite subset of £ and C is a cut of E. Since there are only 
countably many such pairs, Lemma |8~51 implies that there exists an infinite index set Z C {1,2,...} such 
that, for each pair (E, C), 

P E (C):= lim QP(C) 

is well-defined. By Lemma |8~4*1 each Pg is a cutting scheme and there exists a cutting scheme P of £ that 
is consistent with each Pe- 

It remains to verify that P has maximum density r sup with respect to £. Observe that, for any edge 
e G £ and any cutting scheme Q, Q{ e \ is a cutting scheme that chooses among at most two cuts, the set {e} 
itself and (possibly) the empty set, and that the probability that Q{ e } chooses the set {e} is the probability 
that Q chooses some cut that contains e. 
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For any edge e € £ , 



Pr[Pcutse] = Pr[P{ e } chooses the cut {e}] 
= Og}({e}) 

n£i,n-too 1 J 

lim Pr[Q< n ) cuts e] 

n£Z,n— +00 

< lim |e|r sup 

n€l,n- >oo 

Hence, for any edge e, the probability that P cuts e is at most |e|r sup . q 

Lemma 8.7 Lef P be the cutting scheme of £ of maximum density t(£) (shown to exist in Lemma IS!qT ). 
There is a cutting scheme P of 8 that has maximum density t(£) with respect to £. Thus t(£) = t(£). 

Proof: For each point p in the simplex, fix a sequence of rational points (p^^p^, . . .) converging to p. 
For any set F C £ , let denote {{p {n \ q {n) ) : {p, q) G F}. 

Define a sequence of cutting schemes (pW : n = 1, 2, . . .) of £ by p( n )(C) = P(C( n )). That is, P^ 
maps each real point to its nth rational approximation, then cuts using a random cut from P. We claim that 
for each finite E C £ and each cut C of P, the limit 

lim PP{C) = lim P E(n) (C^) 

n— >oc n— >oo 

is well defined. 

To show the claim, it is enough to show that for all e > 0, \P E(n) (C( n )) - P^ (C( m ))| < e for all 
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large enough n and m. Let E' = U E {m ^ and C = U C^ m \ By consistency, 

Pe>{C) < P eM {C^) < P e ,(C')+J2Pe>(D), 

D 

where D ranges over all the cuts of E' other than C that are consistent with C. Each of the cuts that D 
takes on cuts at least one of the segments (p^ , p( m ^ ) where p ranges over the endpoints of the edges in E, 
so the corresponding term in the sum is at most r{£) 

\p(n) _ p M|. Since there are finitely many terms and 
t(£) is finite and p( n > — > p, it follows that the sum on the right is arbitrarily small for large enough n and 
m. Thus \P E ( n) (C^) - P E '(C')\ (and likewise \P E(m) (C^) - P E '{C')\) tends to zero for large m and 
n. Thus so does \P E ( n ){C^ nS} ) — P E ( m ) (C^ m ')\. This proves the claim. 

Define P E (C) := lim^^oo P^T^C). By Lemma 7.4, there is a cutting scheme P of £ that is consistent 
with each P E . The argument used at the end of the proof of Lemma l8^6l shows that, for any edge e € £, 

Pr[P cuts e] = 

< 



Thus the maximum density of P is at most t(P). □ 



lim PH(C) 



lim P f , w (CW) 
lim PrLP cuts e (n) l 

n^oo 

lim r(P)|e (n) | 

n^oo 

t(P) lim |e (n) | 

n— *oo 

r(P)\e\ 
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